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Abstract. Nonlinear ionic screening theory for heterogeneously charged spheres is 
developed in terms of a mode-decomposition of the surface charge. A far-field 
analysis of the resulting electrostatic potential leads to a natural generalization of 
charge renormalization from purely monopolar to dipolar, quadrupolar, etc., including 
'mode-couplings'. Our novel scheme is generally applicable to large classes of surface 
heterogeneities, and is explicitly applied here to Janus spheres with differently charged 
upper and lower hemispheres, revealing strong renormalization effects for all multipoles. 



Screening of heterogeneous surfaces: charge renormalization of Janus particles 2 
1. Introduction 

The past few years have seen an explosion of newly synthesized colloidal (nano-)particles 
that are not spherically symmetric, either by shape (e.g. dumbbells, snowmen, cubes) 
or by surface pattern (patches, stripes) pQ. A broken rotational symmetry also occurs 
when particles are adsorbed to air- water or oil-water interfaces [2], as in for example, 
colloidal monolayers [3j 0], Pickering emulsions [5], or bijels [6]. Moreover, recent 
atomic force microscopy studies have shown that even supposedly homogeneous colloidal 
surfaces can actually be heterogeneous on length scales as large as 100 nm [7], while 
atomic corrugations and facets render any nm-sized particle strictly heterogeneous. An 
important consequence of surface heterogeneity is anisotropy of the mutual effective 
forces, which directly affects the self-assembly process of the (nano) colloids into large- 
scale structures, for instance into ill- understood linear chains [U [7J [8] . 

A fundamental problem is thus to establish relationships between shape and surface 
heterogeneity on the one hand and effective interactions and large-scale self-assembly 
structures on the other [I]. Apart from specific forces (e.g., hydrophobic, van der Waals) 
the effective interactions between dispersed particles often involve a strong generic 
electrostatic component, which is well-described, for homogeneously charged objects, 
by linear-screening theory provided renormalized charges instead of bare charges are 
used [9]. Renormalization of heterogeneously distributed surface charge is an open 
problem, for which we develop a systematic theory in this paper. We go beyond 
recent linear screening treatments [TT] and formulate a new and efficient framework for 
computing nonlinear ionic screening effects of heterogeneously charged spheres dispersed 
in a 1:1 electrolyte. Our theory generalizes Alexander's notion of ion-condensation 
induced charge-renormalization p] to include not only the monopole but also the dipole, 
quadrupole, etc., as well as their nonlinear couplings. These multipole modes can be 
important if one wishes to calculate the electrostatic force between particles, which was 
already shown for clay platelets [10]. Our scheme is versatile and can be applied to 
essentially any type of charge heterogeneity. We focus on applications to Janus spheres 
composed of two differently charged hemispheres [12J. 

We consider an index-matched suspension of N colloidal spheres of radius a in a 
bulk solvent of dielectric constant e and volume V at temperature T. The solvent also 
contains point-like monovalent cations (charge +e) and anions (charge — e) at fugacity 
p s . Here e is the elementary charge. A relatively simple treatment of this many-body 
problem is the cell model [9], in which a single colloid is considered in the center of a 
spherical cell of radius R and volume (47r/3)i? 3 = V/N. We denote the surface charge 
density of this central colloidal particle by ea(9, (p), where 6 and <p are the standard polar 
and azimuthal angle, respectively, with respect to a laboratory frame. Within a mean- 
field approximation, the concentration profiles of the cations and anions can be written 
as Boltzmann distributions p±(r) = p s exp[=)=$(r)], where /csT$(r) / e is the electrostatic 
potential at r = (r,9,(p), with ks the Boltzmann constant and r = |r|. Note that 
$(r) = in the salt reservoir, and that p±(r) = for r < a due to hard-core exclusion. 
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The potential must satisfy the Poisson equation V 2 $(r) = — 4iiX B (p + (r) — p_(r)), where 
we defined the Bjerrum length Xb = e 2 /ekgT '. Combining the Poisson and Boltzmann 
equations gives 

r r 2 I K 2 sinh$(r) r > a, 

where a prime denotes a radial derivative, C? = — [(sin 6)~ 1 do sin Ode + ( srn ^^c^,] the 
angular momentum operator, and = {StxXbPs)^ 1 ^ 2 the screening length. On the 
colloidal surface, r = a, Gauss' law imposes the boundary condition (BC) 

]im$'(r,0,¥>) = lim$'(r,0,^) -47rA B 0-(0,^). (2) 

Electro-neutrality of the cell imposes J dtpdO sin 0$>'(R,0,tp) = at r = R, which is a 
sufficiently stringent BC to close the system of equations in the spherically symmetric 
case. Now, however, an additional BC is to be specified for the angular dependence at 
r = R, depending on the environment of the cell. For now, we assume an environment 
that is characterized by 'isotropic' boundary conditions, 

dg$(R,6,<p)=d ip $(R,9,ip) = 0. (3) 

We will discuss this choice, and its consequences, in section HJ 



2. Theory 

For a given a(9,<p) one can solve ([1]) with BC's for $(r), for example, numerically 
on a discrete (r, 6, <p) grid. The approach we take, however, avoids a cumbersome 
3-dimensional grid in favor of a systematic expansion of the angular dependence in 
spherical harmonics. For notational convenience and illustration purposes we restrict 
attention here to ^-independent cases where the expansion involves only Legendre 
polynomials Pt(x) with x = cos8. 

The first step in this analysis is the decomposition of the colloidal surface charge 
into surface multipoles a> = J^_ 1 dxa(x)Pe(x), such that o~(x) = J2iLo a i^( x )- 
Similarly we decompose $(r, x) = YlT=o^^ r )P^ x )- With $g(r) = (r /a)*$f(a) the 
regular solution to (TjQ) for r G [0, a], the BC's for r G {a, R}, 

$^(a) = -$ e (a) - 4nX B a e (£>0); (4) 

(X 

$o(i?) = and ® t {R) = (£ > 1), (5) 

conveniently decouple for the different £'s. By contrast, the nonlinear sinh term in (TjQ) 
induces 'mode-coupling' between all Legendre components $g(r) - not to be confused 
with mode-couplings in dynamical slowing down. This coupling is obviously unpractical 
for a numerical treatment. 

The second step of our analysis resolves this mode-coupling problem by 'ordering' 
the modes systematically. We introduce a dimensionless 'switching' parameter A, and 
consider the auxiliary distribution a^ A \x) = J^'eLo ^aeP^x), such that A = describes 
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a homogeneous distribution and A = 1 the heterogeneous one of interest. We also 
define the corresponding auxiliary potential & A \r, x) = Xln=o A n <t> n (r, x), where L 
sets the order of the truncation and where the expansion coefficients can themselves 
be expanded as (f> n (r,x) = Y^i=o fnt.{ r )Pi{ x )- The functions f n e(r) are independent 
of A and will be calculated numerically below for the cases of interest n > i (which is 
assumed implicitly from now on). Since the problem is invariant under the simultaneous 
transformation A — > —A and x — > — x one checks that f n i{ r ) — for n + 1 odd, that is, 
we only consider i and n both even or both odd. 

Replacing a(x) by a^ A '{x) and $(r, x) by $^(r, x), inserting the corresponding 
expansions into the BC's, and equating all orders of A yields at r I a and r = R 

£ 

f' nl (a) = -fru(a) - 4:ir\ B o- e 5 n f, (6) 

f nl {R) = (£ = 0) and f n£ (R) = (£ > 1), (7) 

where 5 n e is the Kronecker-delta. When the same replacement and expansion procedure 
is applied to the PB equation (TjQ), one finds upon expanding the argument of the sinh 
term with respect to A a hierarchy of second-order differential equations for f n i(r), 
with a structure that allows for an order-by-order sequential solution. For n = i = 
we obtain for r G [a, R] the spherically symmetric nonlinear PB equation in the cell, 
/oo( r ) + 2 ^°°^ = k 2 sinh foo(r), which we solve explicitly with the BC's given in §6§ and 
j7} on a radial grid. We thus consider /oo( r ) as a known function from now on. For 
n > 1 we obtain 

/» + - (^J^ + cosh/ 00 (r))/ n£ (r) = K 2 S n£ (r), (8) 

where S n g(r) acts as a source term of the form Sn = 0, S20 = \ft\ sinh/oo; 'S'22 = 2S20, 
531 = /u(/20 + 1/22) sinh /00 + 55/11 cosh f 00 , S 33 = | /1 1/22 sinh / 00 + j E f? 1 cosh / 00 , 
and explicit expressions for higher-order terms can be generated straightforwardly. The 
key observation is that S n e only depends on / n /£/'s with n' < n, that is, a hierarchy 
of terms follows spontaneously. Thus ([8]) with the BC's ([6]) and (J7J) can be solved for 
n = £ = 1, which in turn determines S2o( r ) an d 822(1") such that /2o( r ) and can 
be solved, etc. The nonlinear mode-coupling, represented explicitly by cosh/ o(r) and 
S n i(r), renders the linear equation ([8]) highly nontrivial, yet numerical solution on a 
radial grid r e [a, R] is straightforward. With f n t{f) determined for L > n > £ > 
for some cut-off L, we can set A = 1 to explicitly construct the potential of interest 
*(r,x) = J2toMr)Pt(x) with $,(r) = EtuMr). 

3. Results 

The theory developed so far is directly applicable to any uniaxial charge distribution, 
while generalizations to azimuthal dependencies and nonspherical shapes are feasible. 
In this paper we illustrate our scheme for the prototype heterogeneous charge distribtion 
of Janus spheres characterized by surface charge-densities o~n and as on the northern 
(x > 0) and southern (x < 0) hemisphere, respectively [T2l [13]. The non- vanishing 
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Figure 1. (a) Angular dependence of the potential (at several distances r) and the 
surface charge (inset) of a Janus sphere of radius a, for na — 1, ZqXb/cl = 22.5, 
R/a = 3, and L — 13, for charge heterogeneities (see text) g = (arrows), g = 0.5 (dot- 
dashed), and g = 0.8 (solid curves). The dotted curves for r/a = 1.5 and g S {0.5, 0.8} 
stem from a Yukawa-segment model (see text), (b) Radial dependence of the monopole 
(£ — 0), dipole (I — 1), and quadrupole (I = 2) potentials $^(r) (solid lines), and their 
far- field linear-screening approximations (dashed curves), for the parameters of (a) and 
g = 0.5. The angles defined at r = a relate to the bare and renormalized modes of the 
surface charge (see text). 



modes are thus cr = (a^ + 0s)/2 and o~i = ga (2£ + 1)(— 1) 2 for £ odd. Here we 

defined the dimensionless heterogeneity parameter 

O-N-O-s , , 

9 = rr M„ ' ^ 

a N + a s 

which together with the total charge Zq = 4ira 2 ao fully characterizes the distribution. 
Below we set na — 1, R/a — 3, and L — 13 throughout unless stated otherwise, 
and we identify ZqXb/o as the only other relevant dimensionless combination. For 
Z \b/o = 22.5 and g G {0,0.5,0.8} figure 1(a) shows the ^-dependence of a (inset) 
and $ for several r, revealing isotropy (as expected) for a homogeneous surface charge 
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(g = 0, arrows) and strong anisotropy for the heterogeneous cases g = 0.5 (cr/v — 3(75, 
dot-dashed) and g = 0.8 (o~n = 9o~s, solid lines); the (small) oscillations with 9 at 
r = a are numerical artefacts due to the truncation at L = 13. The ^-dependence of 
the potential weakens, as expected, for increasing distances r. Figure 1(b) shows $^(r) 
for £ = 0, 1, 2 and g = 0.5. Interestingly, the modes with £=1,2 have a non-vanishing 
electric field in the interior of the particle. The overall magnitude and spatial variation 
of $ in figure 1 show the need for nonlinear screening theory. Nevertheless, in analogy 
to the spherically symmetric case [9] one can describe the far-field potential (r ~ R) 
and hence the colloidal interactions in terms of linear screening theory (dashed curves 
in figure 1(b)) with a renormalized surface charge distribution o~*(x) = Ylt a e^( x ) that 
we will calculate below. 

In the far-field r ~ R we treat the deviation of $(r, x) from its angular average 
$ (i?) = Xo a t r — R as a small expansion parameter, such that (0Q) for r > a can be 
linearized as V 2 $(r, x) ~ K 2 [tanhx + (^K r > x ) ~ Xo)] with R 2 = /t 2 coshx - The uniaxial 
solutions to this linear PB (LPB) equation read $(r, x) ~ xo — tanhx + X^o[ a ^(^ r ) + 
&^(/«r)]P^(:r) where it and fc^ are modified spherical Bessel functions. The coefficients 
at and bt are integration constants that we fix by matching the LPB-solution at r = R, 
for each £, to <&e(R) = an d ^(-R) = °f the nonlinear problem. This leads for every 
£ to the linear two by two problem 

Xe = (Xo - tanhxo)^o + aik{RR) + b e ki(RR); 

X' £ = R\ati' t (RR) + b e k' e (RR)^j, (10) 
which results in explicit expressions for at and bt given by 

a _ vt (k e+ i(RR) - jjjk e (RR)) + x'thi^R) _ 
ii(RR)ke + i(RR) + ie+i(RR)ke(RR) ' 

b ^ = vi (U+i(kR) + -j^iiRR)) - x'M^R) > 12 v 
ie(RR)ki + i(RR) + i£ + i(RR)ke(RR) 

where vt = Xe~ (Xo ~ tanhx o )^ . The dashed curves in figure 1(b) are the result of such 
a far- field fit. With at and bg explicitly known, one can extrapolate the LPB solution 
to r = a to yield, with (j3J) and standard Bessel function relations, the renormalized 
multipoles 

<y\ = - 4~^^( a ^+ l( ^ a ) - b e k e+1 (Ra)y (13) 
This expression is the multipole generalization of the well-known charge renormalization 

A first illustration of multipole renormalization is shown by the angles in figure 
1(b), which represent slope discontinuities at r = a which are proportional to at (non- 
linear theory) and a\ (far- field fit). For £ = 0, 1 we see a\ < at, which means that the 
effective charge and dipole are renormalized downward. Interestingly, however, for I = 2 
the inset in the lowest panel reveals upward renormalization since o~\ ^ while a2 = 0, 
that is, the Janus particle has a mode-coupling induced far-field quadrupole signature. 
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Figure 2. Renormalized (scaled) monopole Zq for several g, L (see text) of a Janus 
sphere as a function of the total charge ZqXb /a. The graph includes data for cell radii 
R/a = 1.5 and R/a = 3.0. The dot-dashed curve denotes Zq = Zq. 
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Figure 3. Renormalized (scaled) dipole Z\ for several g, L (see text) of a Janus sphere 
as a function of the total charge ZqXb/cl- The upper and the lower graph represent 
cell sizes R/a — 3.0 and R/a = 1.5 respectively. The thick dotted lines denote 
the (scaled) difference of the renormalized northern and southern charge presumed 
distributed homogeneously (see text) for g — 0.8. The dot-dashed curves in both 
graphs denote Z\ = Z%. 



Multipole renormalization is quantified further in figures 2-4, where (scaled) renormal- 
ized multipoles Z\ = 4na 2 a^ are shown for the monopole £ = in figure 2, the dipole 
I = 1 in figure 3 (scaled with g), and the higher-order multipoles i = 1, ■ • • , 6 in figure 
4. All multipoles are shown as a function of Z Xb/o. Furthermore, in figure 2 and 3 we 
chose several truncation levels L and heterogeneities g, and also picked two cell radii R. 
The dot-dashed curves denote the linear limit Z\ = Zg. Figure 2 and 3 show that all 
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Figure 4. Rcnormalized (scaled) multipoles Z* t of a Janus sphere as a function of the 
total charge Z^Xs/a. The parameters here are g = 0.4 and L = 13. 

curves for L = 13 superimpose on those of L = 9, 11 for all Zq and g < 0.8, indicative of 
excellent convergence in this parameter regime; for g = 1 the convergence deteriorates 
for ZqXb/o > 15. The R dependence in figure 2 and 3 shows the strongest renormal- 
ization in the largest cell, not unlike the homogeneous- charge case [9]. Interestingly, in 
the nonlinear regime Z X B /a > 10 figure 2 shows a mode-coupling induced reduction 
of Zq by 10 's of percents when g increases from to 1. In other words, in contrast to 
the more usual 'linear' electrostatics we now have a far-field monopole potential that 
is not only determined by the net charge but in fact also by its heterogeneity. This is 
a key finding, relevant for understanding patchy-particle interactions. The mode cou- 
pling has an even stronger effect on renormalization of Z*, for which figure 3(a) and 
(b) show a pronounced maximum in between the \ow-Zq linear screening regime and 
the high-Z regime in which Z\ becomes even vanishingly small for all g's considered. 
Figure 4 shows, for g = 0.4, that in fact all Z\ with i ^ vanish in the limit of large 
Z , while they all show an intermediate regime with finite values even for i = 2, 4, 6 
for which o~e = 0. The underlying physics for non-oppositely charged hemispheres with 
o~n > o~s > (that is, with < g < 1) is that both a^ and as renormalize, if both 
are high enough, to the same saturated value, giving rise to a pure far-field monopole 
without multipoles. 

The idea might emerge that both hemispheres renormalize their charge independent 
of each other, that is, some of the results could suggest that the renormalized charge 
density on the colloidal surface is a function which depends only locally on the bare 
charge density. If that were the case, it would suffice to calculate the renormalized 
surface charge- density for a^ and as, as if both were the charge density of a monopole. 
The cell radius is to be kept unchanged. The thick dotted lines in figure 3(a) and (b) 
denote the resulting dipole charge for g = 0.8, which is calculated with the obtained 
(renormalized) a* N and a* s by \{a* N — a* s ). The correspondence with the solid line is at 
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best reasonable but not perfect. Also the effective monopole is not accurately predicted. 
This can can be seen from figure 5, where we investigate the cell-size dependence of the 
the effective monopole for g = 0.4, 0.8 and for Ka = 1, 3. We included data from the full 
theory (depicted by symbols) and the predicted values by 2na 2 (a* N + a* s ) as a solid line. 
One can see that the monopole charge is underestimated for a wide range of cell radii, 
especially for larger cells. Nevertheless, there is qualitative agreement on the increase of 
the renormalized charge with higher values for Ka. The difference with the full theory 
is expected to be the largest for very heterogeneously charged particles. Indeed, we see 
the largest discrepancy in figure 5(a) and (b) for g = 0.8, with deviations up to 20-25% 
between the results of the present theory for ctq and those of the simple approximation 
(cr^r + o~g)/2 discussed above. Apparently, the interactions between the hemispheres do 
play a role, which in fact can can also be concluded from the induced even multipoles 
in figure 4. Further research might give more insight into the characteristics of these 
interactions. 
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Figure 5. Renormalized monopole charge Zq for several g of a Janus sphere as a 
function of the cell size R, both for ZoXs/a = 22.5, in (a) for Ka — 1 and in (b) for 
na = 3. The thick solid line corresponds to a pure bare monopole, yielding g = 0. 
The diamonds (g = 0.4) and circles (g — 0.8) show the data from the full theory, 
using L — 13. The thinner solid lines denote the average of the renormalized northern 
and southern charge presumed distributed homogeneously (see text) for g = 0.4 and 
g = 0.8. 



4. Discussion and Conclusion 

The newly emerging general picture is that nonlinear ionic screening of heterogeneously 
distributed surface charges strongly affects the far-field symmetry of the potential, and 
hence also the symmetry of the effective interactions and the self-assembling structures. 
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The systematics of the present screening theory could be a firm basis to further study 
these intricate features of heterogeneously charged particles. The category of particles 
which are described by g > 1, carrying positive and negative charges on the two 
hemispheres, are particularly interesting. Because the present method is not very 
efficient at high values of g, we have developed an alternative method to treat these 
mainly dipolar particles within a mode expansion as well. This will be presented in 
further work. 

It is tempting to model Janus spheres by a Yukawa-segment model p3] in which 
every surface element dS contributes a* ^Xbs' 1 exp(—Ks)dS to the (dimensionless) 
potential at a distance s. Here the renormalized charge densities on both hemispheres 
is obtained from the renormalized monopole and dipole charge density via 



The dotted curves in figure 1(a) show that agreement with our full calculations 
is reasonable though not quantitative; the Yukawa model ignores the ionic hard- 
core exclusion in the interior of the particle. Therefore, effectively it describes the 
(dimensionless) potential of a charge configuration in which oppositely charged ions 
were able to approach the heterogeneously charged surface from two sides, such that 
this potential is more suppressed compared to the full theory. 

Being a point of discussion, we return to the choice of the boundary conditions 
on the cell's surface ([3]), which we called 'isotropic' BCs. This denomination follows 
from the fact that, considering two randomly oriented neighboring cells, the cell-surface 
potentials of two cells should match on the spot where they touch, giving rise to a 
constant cell-surface potential. Nevertheless the choice of BCs is not unique. We 
can also supply the system with 'nematic' boundary conditions, corresponding to the 
situation that all cells are perfectly aligned such that cells only touch on opposite 
spots. On these spots the electrostatic potential should match, and we can also demand 
continuity of the electric field. The BCs then become 



In fact one can even interpolate between 'isotropic' and 'nematic' BCs by introducing 
an orientation distribution function [14] . In this article, we assume a system in which 
the cell boundary is best described by isotropic BCs, given by ([3]). However, for the 
parameters used in this article, it turns out that this particular choice for the BCs did 
not noticeably affect the values of the renormalized charges. We do not see a significant 
change by turning to nematic BCs ffl5|) . This insensitivity to the choice of BCs is due 
to the fact that the nonlinear behaviour is an effect which takes place close to the col- 
loidal surface, where these BCs have the least influence on the electrostatic potential. 
Furthermore, the only multipoles which are directly affected by the particular choice of 
BCs are the nonzero (£ > 0) even multipoles, which are small for Janus particles. We 




(14) 



$(R,8,<p) = $(R,ir-9,ir + <p), 
$'(#, 9,(p)= - $'(R, it-6,it + v). 



(15) 
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therefore think that the obtained values for the renormalized multipoles can be applied 
in a model to describe the behaviour of a many-body system within linear theory, no 
matter what the orientations of the surrounding colloids are. Since the monopole and 
multipole potentials decay for large r equally fast as exp(— Rr)/r [15], where k -1 is the 
decay length, the renormalized multipole charges are expected to contribute in dense as 
well as dilute systems. 

In summary, we have developed a systematic framework for nonlinear ionic 
screening of heterogeneously charged spheres. The scheme allows for an explicit far- 
field analysis that generalizes charge renormalization from the well-studied homogeneous 
case (pure monopole) [9] to the heterogeneous case (dipoles, quadrupoles, etc. and 
their nonlinear couplings). Application to charged Janus spheres shows (i) a 40% 
reduction of the effective monopole for g = 1 (charged and uncharged hemisphere) 
compared to g — (homogeneously charged sphere), (ii) a mode-coupling induced 
far-field effective quadrupole component without an actual surface quadrupole, (iii) 
a pure far-field monopole with vanishing higher-order multipoles in the saturated 
high-charge limit, and (iv) no quantitative agreement with a simple Yukawa-segment 
model based on renormalized multipoles. Our study opens the way to systematic 
microscopic calculations of effective electrostatic interactions between Janus (and other 
patchy) particles. In addition, our analyses also reveal non-vanishing electric fields 
inside heterogeneously charged particles, which could couple to interior dipoles and 
affect (anisotropic) mutual Van der Waals forces. Given that the presently introduced 
expansion technique can be generalized to other geometries (for example, patterned 
planar surfaces or ellipsoidal patchy colloids), our technique and findings are directly 
relevant for gaining microscopic understanding of effective interactions and ultimately 
phase behaviour of a large class of dispersions of patchy or patterned nanoparticles, 
colloids, or proteins pQ. 
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